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We present theoretical results for the backaction force noise and damping of a mechanical oscillator 
whose position is measured by a mesoscopic conductor. Our scattering approach is applicable to a 
wide class of systems; in particular, it may be used to describe point contact position detectors far 
from the weak tunneling limit. We find that the backaction depends not only on the mechanical 
modulation of transmission probabilities but also on the modulation of scattering phases, even in 
the absence of a magnetic field. We illustrate our general approach with several simple examples, 
and use it to calculate the backaction for a movable, Au atomic point contact modeled by ah initio 
density functional theory. 
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Quantum mechanics requires that any detector used to 
measure an object's position unavoidably exerts a back- 
action force, imposing a fundamental limit on continu- 
ous position detection [1, 2]. Recent experiments with 
nanoelectromechanical systems (NEMS) have come re- 
markably close to realizing this limit by using quantum 
electronic conductors as position detectors of nanome- 
chanical oscillators [3-5]. In these systems, position de- 
tection is achieved using the influence of the mechanics 
on the current through the conductor; thus, it is natural 
to associate backaction with the position sensitivity of 
the electron transmission probability. This is indeed the 
picture that emerges from theoretical studies in the limit 
of weak tunneling [6-8] ; however, several recent experi- 
ments are far from this limit [5, 9-11], and it is not clear 
that the weak tunneling results apply. 

In this paper, we study the backaction of a mesoscopic 
position detector using a general noninteracting scatter- 
ing approach that is not limited to the weak tunneling 
limit. Scattering theory has been used extensively to 
study various aspects of mesocopic conductors, and we 
adapt it here to study the backaction heating and damp- 
ing of a mechanical oscillator coupled to a conductor. 
Surprisingly, we find that backaction arises not only from 
transmission probabilities, but also from the position sen- 
sitivity of scattering phases, and present several simple 
but illustrative examples where the phases play a pro- 
nounced role. We emphasize that these phases may be 
important despite intact time reversal symmetry, which 
we assume throughout, unlike Aharanov-Bohm phases 
due to a magnetic field [8] . Finally, we apply our general 
results to calculate the backaction from an atomic point 
contact (APC) between Au electrodes, using a scattering 
matrix obtained from density functional theory (DFT). 

Our approach significantly extends the seminal work of 
Yurke and Kochanski, who first considered force noise in 
a tunnel junction using a scattering approach [12]. Un- 
like their work, which is limited to particular scatter- 
ing potentials, we rely only on general properties of the 
scattering matrix. As a result, we can describe a wide 



class of systems including arbitrary scattering potentials, 
various forms of electromechanical coupling, and multi- 
channel scatterers. Moreover, we calculate not only the 
backaction force noise but also the backaction damping, 
which is important in experiments (e.g. it is the basis of 
backaction-cooling [4] ) and up to now has not been dealt 
with in the scattering approach. 

Scattering approach. — We consider a two-terminal de- 
vice consisting of a coherent scattering region coupled to 
left and right leads, each of which supports N transverse 
modes. Electrons are scattered by a potential U(r e ,x), 
where r e = (x e ,y e ,z e ) is the electron position, and the 
potential depends on the position x of a mechanical os- 
cillator. Incoming and outgoing waves are related by 
the scattering matrix s(x), which depends on x through 
U(r e ,x). We will show that a knowledge of s and ds/dx 
is sufficient to calculate the backaction. 

For the usual experimental regime of weak electrome- 
chanical coupling, the change in the electronic potential 
due to small changes in x are generally linear and may be 
written Hint — —xF, where the force on the oscillator is 
F = — J dr e p(r e )dU(r e ,x)/dx, and p(f e ) is the electron 
density operator. By relating small and slow changes in 
the potential, U(r e ,x), to the parametric derivative of 
the scattering matrix [13-15], we can express F in the 
scattering state basis as 




where a a (e) destroys a scattering state of energy e inci- 
dent in lead a, and 

We require only the diagonal-in-cncrgy part of W since 
we focus on the zero frequency noise properties of F, 
sufficient for the experimentally relevant case when the 
oscillator period is much longer than timescales in the 
conductor. Derivatives of the scattering matrix similar to 
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Eq. (2) are familiar from studies of charge noise [16] and 
parametric pumping [17]. Here we use the parametric 
derivative with respect to x to calculate the backaction on 
the oscillator directly in terms of the scattering matrix, 
without the need for a detailed knowledge of U (f e , x) and 
p(f e ) in the scattering region. In the following we work 
to lowest order in H ln t, valid for weak coupling. 

Fluctuations of the backaction force cause momentum 
diffusion and heating of the oscillator. Heating is de- 
termined by the classical, frequency-symmetric part of 
the backaction force noise, Sf[u] — (Sf[u] + Sf [—<*>]) /2, 
where the quantum noise spectral density is Sf[w] = 
J dte lu>t (F(t)F(0)) and averages are taken with respect to 
the uncoupled conductor [18]. These averages are easily 
taken using Eq. (1), and the backaction heating is di- 
rectly determined by W. The zero frequency force noise 
is (k B = 1,S F = S f [0}) 



S f = 2-kTi 



ap 



debr{W aP W Pa }f a (l-fp), (3) 



where the trace is over transverse modes, assumed to 
be the same in both leads, and the matrixes W a p are the 
NxN blocks of W in Eq. (2) , which may be e-dependent. 
The Fermi functions are f a = (1 + e f £ "' lt, '/ Tcl )~ 1 , where 
fj, a is the chemical potential in lead a and T e \ is the elec- 
tronic temperature . 

In addition to heating, the oscillator also experiences 
backaction damping as a result of energy exchange with 
the conductor. The damping rate is given by the quan- 
tum, asymmetric-in-frequency part of the force noise, 
7[w] = (SVM - S f \-u\) /2MTioj, where M is the os- 
cillator mass. Taking the u — > limit, we find 
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(4) 



By considering the ratio of Sf[w] to 7[w], one can 
associate a frequency-dependent effective temperature, 
r c ff[w], with the backaction; this amounts to using the 
standard fluctuation-dissipation relation to define the ef- 
fective temperature at each frequency from the system's 
force noise and damping [2, 6, 18]. TcffM characterizes 
the conductor as an effective thermal environment. In 
the uj — > limit, the relation is simply T e g = Sp/2Mj. 
If backaction dominates over intrinsic sources of dissipa- 
tion, T ff corresponds to the physical temperature of the 
oscillator. 

Single channel. — We first consider the case of single- 
channel leads. For simplicity, we also focus on the 
limit of small applied bias, ignoring the possible energy- 
dependence of s. We assume time reversal symmetry 
(i.e. no magnetic field), but allow for broken left-right 
inversion symmetry. In this case the scattering matrix 
may be parametrized as 

Vne ie iVT \ 



s(e, x) = e 



(5) 



where T (TZ = 1 — T) is the transmission (reflection) 
probability, <f> is the overall scattering phase, and 
parametrizes broken inversion symmetry, i.e. 9 = for 
an inversion-symmetric conductor. In general, all of the 
scattering parameters depend on x through the potential, 
U(f e ,x). 

Inserting Eq. (5) into Eq. (3) we obtain the sym- 
metrized force noise for a single channel, 



S F = ^ &r ' d feV 
2tt 4TZT 



(1 + TIT A g ) coth 



(6) 
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where V is the bias, and the phase terms enter as 



A c = WIT 



( dC/dx 
[dT/dxJ ' 



(7) 



for £ = <j>, 6. In the limit eV 3> T c \, the first term (in- 
dependent of and Ag) in Eq. (6) represents the ex- 
pected, quantum-limited backaction of our position de- 
tector: it is simply the sensitivity of a position measure- 
ment by monitoring the current, and reflects the fact 
that a stronger measurement leads to increased backac- 
tion. This term scales as the square of the measurement 
gain, xif « dT/dx, and inversely with the shot noise 
in the current, 5/ = e 2 V7ZT /2irh; in the limit T < 1 
it reproduces the well-known result obtained from a tun- 
nel Hamiltonian calculation [6, 7]. The second term in 
Eq. (6) is independent of dT/dx and is thus not di- 
rectly related to a measurement of the current. Instead, 
it results from the oscillator's modulation of the phase 
0. This phase contribution to Sf is proportional to 1ZT 
and thus vanishes when T< 1. The remaining two terms 
(oc T e \/eV) are also independent of dT/dx and represent 
additional thermal noise at finite T e \/eV. 
The damping for a single channel is 



7 = 



h (dT/dx) 2 
2nM ATIT 



(l + A + ^A e ). 



(8) 



In the small bias limit, 7 is strictly positive and indepen- 
dent of T c \. Similar to Sf, the first term in Eq. (8) is the 
backaction associated with a measurement of the current 
and reduces to the tunnel Hamiltonian result in the limit 
Tel More interestingly, the second and third terms 
correspond to corrections due to scattering phases; un- 
like Sf, these phase contributions to 7 are present even 
for a symmetric detector and, as we will see, do not nec- 
essarily vanish in the weak tunneling limit. The overall 
phase (j> is directly connected to the density of states in 
the scattering region via the Friedel sum rule [19]. An 
x-dependent <fi implies that the mechanical oscillator can 
change the scattering- induced electronic density of states; 
this means that the total electronic free energy becomes 
x-dependent, resulting in a force whose quantum noise 
contributes to damping. 



Equations (6) and (8) show that the backaction prop- 
erties of a general conductor cannot simply be extrap- 
olated from the weak tunneling limit; scattering phases 
play a role in both the heating and damping of the oscil- 
lator. Further, the phases can have a dramatic influence 
on the effective backaction temperature T e g of the detec- 
tor. For a single channel, using Eqs. (6) and (8), in the 
limit T e i < eV we find 



(a) (c) T L (x) erf T R (e) r L (x) ta Mx) 



eff 



eV f 1 + KTAg 
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If the mechanical motion does not modify the scattering 
phases, then we simply obtain the tunnel Hamiltonian 
result [6, 7], T c s = eV/2, independent of T. However, 
in the more general case including the backaction from 
scattering phases, T e s is not solely determined by the 
voltage. The phase corrections always decrease the effec- 
tive temperature; they arise from the diagonal elements 
of W, which correspond to transitions between scatter- 
ing states in the same lead. At T c \ — such transitions 
can only occur if an electron absorbs energy, because the 
scattering states in each lead are filled up to the Fermi 
level. Thus, phase corrections lead to increased absorp- 
tion of energy from the oscillator, lowering T e ff. Including 
a non-zero lead temperature T e \, one finds that T e g can 
be lowered to a minimum value of T c \; as T c \ <C eV, this 
could still be quite useful. 

Square potential barrier. — To demonstrate that back- 
action from scattering phases plays a role even in the 
simplest scattering model, we calculate the backaction 
for a one dimensional symmetric square barrier poten- 
tial whose width depends on the oscillator position. The 
force noise for this model was first considered in Ref. 12; 
our general method further provides 7 and T e g and al- 
lows us to identify the role of scattering phases. Incoming 
electrons of wavevector k and energy e are scattered in 
one dimension by a square potential barrier of height Uq 
and width w = L + x. The inverse decay length of the 
wavefunction under the barrier is n = y // 2m e (Uo — E)/h. 
It is straightforward to find T and <fi as functions of Uq 
and w, and 9 = due to inversion symmetry. We obtain 



2^2 



(fc 2 -K 2 )T, 



(10) 



and the phase terms become important when T ~ 1. For 
a high but narrow barrier {Uq ^> e, kL <C 1), we find (via 
Eq. (9)) that T e g may be reduced by up to a factor two 
compared to the tunnel Hamiltonian result of eV/2. For 
a low barrier (Uq C e), we find T c g — > T c \. 

Resonant level model. — We now apply our general re- 
sults to a prototypical resonant level model (RLM), 
where a single electronic level of energy e d is connected 
to the left (right) lead via tunneling rate Tl (Fr). If 
td depends on the position of a mechanical oscillator (see 
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FIG. 1: Schematic setups and backaction damping for RLM 
with dispersive (a,b), dissipative (c,d) and shuttle (e,f) me- 
chanical coupling. In all plots, the full backaction damping 
(red solid) and the damping without phase corrections (blue 
dashed) are shown, with T (black dotted) also shown for ref- 
erence. We took r L (0) = T R (0). 



Fig. la), one has the electromechanical analog of a disper- 
sively coupled optomechanical system [20] , and a simple 
model of quantum-dot-based NEMS studied in recent ex- 
periments [21, 22]. Beginning from the scattering matrix 
for the RLM, s a/3 = S af3 - ihy/r a Tp/ (e - t d + ihT/2), 
where T = Tl+Tr, and assuming linear coupling, we ob- 
tain 7 cx {T/T L T R ) 2 T 2 [23]. We also find A = l + "KA e , 
independent of the tunneling rates and the detuning of 
the incident electron energy e from e^. Comparing with 
Eq. (8), this implies that the phases play a crucial role: 
the x- dependence of the overall scattering phase <f> always 
accounts for half of the damping, as seen in Fig. lb. Fur- 
ther, in the limit of asymmetric tunneling rates we find 
7?.Ag 3> 1, and the damping is almost entirely due to the 
combined contributions from (f> and 9. Also striking is the 
cotunneling limit, where the detuning is large compared 
to the level broadening, i.e. \e — ed\ 3> KT. In this limit 
tunneling is suppressed, T <C 1, and the level charge only 
fluctuates virtually; as a result, one might expect that the 
system is equivalent to a single junction in the weak tun- 
neling limit, and that 7 should be given by the tunnel 
Hamiltonian result, i.e. the first term in Eq. (8). How- 
ever, this is not the case: due to phase corrections, the 
damping is twice the tunnel Hamiltonian result (see inset 
of Fig. lb). This shows that the phases can play a role 
even when T is small. Note that backaction in this model 
was recently studied theoretically using a path integral 
approach [24, 25], although backaction due to phases was 
not discussed. 

Our general theory also allows us to consider varia- 
tions of the above RLM where the mechanical position 
modulates the tunneling rates T L and T R . This is the 
electromechanical analog of a dissipatively coupled op- 
tomechanical system [26], and could be achieved experi- 
mentally using a quantum dot coupled to two leads via 
tunnel junctions, with the tunneling rates modified by an 
on-board [5] or off-board [10] mechanical oscillator. First, 
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L (A) L (A) 

FIG. 2: (a) APC in a 5 x 5 atom Au (100) quantum wire, 
(b) Conductance through the APC versus gap size (for small 
L, several channels contribute to transport and G exceeds 
e 2 /h [23]). (c) Effective backaction temperature versus gap 
size. For small L, T e « < eV/2 due to phase corrections. The 
Au-Au interatomic distance in bulk gold is 2.87 A. 



transmission properties are often studied using DFT, an 
important feature of our calculation is our explicit use of 
the scattering phases obtained from an atomistic calcu- 
lation of a quantum electronic device. 

Conclusions. — We have presented a scattering ap- 
proach to backaction in NEMS and demonstrated the 
importance of backaction from scattering phases. This 
work is particularly relevant to NEMS based on quan- 
tum or atomic point contacts which are often far from 
the weak tunneling limit. Our results may also be easily 
extended to describe strong electromechanical coupling 
in the low oscillator frequency limit, by making an adi- 
abatic approximation such that the noise spectra of F 
effectively become ^-dependent [18]. 

This work was supported by NSERC, FQRNT and CI- 
FAR. 



we consider a setup where only the left tunneling rate is 
^-dependent (see Fig. lc,d). In this case, interference 
between resonant charge fluctuations (on the level) and 
non-resonant charge fluctuations (in the leads), result in 
a Fano lineshape and suppression of 7 at zero detuning 
[23], similar to the optomechanical case [26]. Second, 
we consider mechanical coupling to both tunneling rates 
with opposite sign, corresponding to a quantum shuttle 
(see Fig. le,f). Here we find 7 cx (T /T L Y R ) 2 T; moreover, 
all of the damping is due to the scattering phases, since 
(for Fi = Tr) the transmission has no linear dependence 
on x [23]. 

Atomistic model. — While the above examples show 
that phases contribute to backaction in simple model po- 
tentials, our approach allows us to investigate phase con- 
tributions in fully atomistic calculations of mesoscopic 
conductors. We demonstrate this by applying our theory 
to an APC using the scattering matrix obtained from 
DFT [27]. We model the APC as a single-atom constric- 
tion in a 5x5 atom Au quantum wire (see Fig. 2), and 
take x to modify the gap size of the APC to L + x. This 
geometry is motivated by recent experimental setups us- 
ing an APC [5] or scanning tunneling microscope (STM) 
[9] with one mechanically compliant electrode. We ap- 
proximate the surface electrode of experiments by the 
flat 5x5 edge of wire on the right; this is justified since 
the transport properties of the APC are expected to be 
dominated by the few atoms closest to the tip. We find 
11 scattering channels contributing to transport, consis- 
tent with recent ab initio studies of similar Au wires [28] . 
After obtaining s and ds/dx [23], we calculate the back- 
action using Eqs. (3) and (4), assuming T cl <C eV. We 
find that phase corrections are important when the APC 
transmission deviates from the weak tunneling limit; it 
leads to a significant reduction in T e g from the tunnel 
Hamiltonian result of eV/2, as seen in Fig. 2. While 
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ERRATUM ADDED 

In Ref. [1] , Eq. (4) for the backaction damping is valid 
only when the scattering matrix s(e, x) is independent of 
energy over the range of energies contributing to scat- 
tering; we erroneously stated that this expression was 



correct for general, energy-dependent scattering. Using 
the definition of the W matrix given in Eq. (1) of Ref. [1], 
one can easily show that the full expression for the low- 
frequency, linear-response backaction damping (valid for 
arbitrary energy-dependent scattering) is 



i = S / de tr £ w °* w e° (- d -k) + 2 to - /*) 



We stress that this correction has no impact on the subse- 
quent results of Ref. [1], as we exclusively considered the 
limit of small drain-source voltages in which the energy 
dependence of the scattering matrix plays no role. 

It is worth briefly outlining the different origins of 
the terms in Eq. (El). The first term, proportional to 
W a pWp a evaluated at energy e, corresponds to the damp- 
ing discussed in Ref. [1]. This term arises from the in- 
crease in the number of scattering transitions contribut- 
ing to the force noise when an electron in the conduc- 
tor absorbs energy u from the mechanical degree of free- 
dom (and conversely, the decrease if an electron emits 
energy to the mechanics). Because of this asymmetry, 
the conductor favors absorption of energy and this part 
of the damping is always positive. In contrast, the term 
in Eq. (El) proportional to /l — Jr constitutes a nonequi- 
librium contribution originating from the energy depen- 
dence of the matrix elements of the force operator be- 
tween scattering states at different energies. This in- 
trinsic energy dependence of the matrix elements may 
favor either absorption or emission of energy resulting in 
positive or negative damping. This contribution was re- 
cently presented in Ref. [2] for a general model where the 
electronic system is a multilevel, noninteracting quantum 




r 



dot. 

It is interesting to note that, following the lines of 
Ref. [3], it is not possible to relate the nonequilibrium 
term in 7 to the "frozen" scattering matrix s(e, x); one 
also needs knowledge of the form of the wavefunctions in 
the scattering region. Nonetheless, one can easily iden- 
tify general classes of systems where it plays no role. We 
find that the nonequilibrium term always vanishes for an 
inversion-symmetric scattering potential, as well as for 
the single-resonant-lcvcl models considered in Ref. [1]. 
The non-equilibrium damping terms thus play no role 
(even at finite bias voltage) for the examples considered 
in Ref. [1]. 

We thank F. von Oppcn for drawing our attention to 
the incorrect presentation of Eq. (4) in our original paper. 
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RESONANT LEVEL MODEL 

As given in the main text, the scattering matrix for 
the RLM (see, for example, [1]) is 



s af )(e,x) = 8 a p - 



ir/2' 



(SI) 



where e d is the level energy, T L (T R ) is the tunnel- 
ing rate to the left (right) lead and r = Tl + T^. 
Here we set Ti = 1. The transmission probability is 
T = T L Y R / [(e - e d )) 2 + T 2 /4] . In Fig. 1 of the main 
text we sketch and plot the damping for three forms 
of mechanical coupling. First, dispersive coupling corre- 
sponds to an .x-dcpendent level energy and we replace e d 
with id(x). We assume linear coupling, id(x) = €d + Ax, 
and obtain 



A 2 



AttM \( e - 



[(e-e d ) 2 +r 2 /4] 



2 • 



(S2) 



and we see that 7 cx {Y/TtXr) / as stated in the 
main text. Second, dissipative coupling corresponds to 
^-dependent tunneling rates. We take the left rate to de- 
pend linearly on x and replace Tl with Tl(x) = Tl+Ax. 
This yields 
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a 2 (r + r L )(6- £d ) 2 + r fl r 2 /4 

4ttM 2r L [(e-e d ) 2 +r 2 /4] 2 



(S3) 



This result can be verified by a direct calculation start- 
ing from the Hamiltonian for a single resonant level, and 
using the same Tl(x) [2]. Finally, we model a quantum 
shuttle with ^-dependent left and right tunneling rates, 
f L (x) = T L + Ax and f ^(a;) = T R - Ax. We find 



7 = 



_^ r 2 

A^MT L T R [{e-e d ) 2 + T 2 /AY 



(S4) 



and we see that 7 cx (r/r£r^) 2 T as stated in the main 
text. 



ATOMIC POINT CONTACT MODEL 
Atomic structure and ab initio methods 

The APC consists of a tapered 5x5 atom Au(100) 
quantum wire opposing an identical flat-surfaced 



nanowire (see Fig. 2 in the main text). The atomic posi- 
tions of the three Au layers forming the tip (on the left) 
and the first Au layer of the flat surface (on the right) 
were independently relaxed to forces < 0.01 eV/A us- 
ing the VASP density functional theory (DFT) software 
package [3, 4]. All non-relaxed atoms were fixed to their 
DFT optimized bulk positions corresponding to a lattice 
constant of 4.06 A for FCC gold. 

From the optimized APC structure, we performed the 
first principles transport calculations using the MatD- 
Cal device simulator, based on a combination of DFT 
and Keldysh nonequilibrium Green's functions (NEGF) 
[5, 6]. MatDCal uses a linear combination of atomic or- 
bitals as a basis, where two s-orbitals, one p-orbital and 
two d-orbitals were employed for each Au atom. The 
exchange and correlation energies were described within 
the local density approximation [7] , while the nuclear and 
core electrostatic potentials were modelled using norm- 
conserving non-local pseudopotentials [8]. Note that the 
size of the supercell box, in the directions perpendicular 
to the transport direction, was chosen large enough to en- 
sure no interactions with the neighboring supercells. Af- 
ter calculating the self-consistent DFT Hamiltonian, we 
solve for the scattering matrix which provides the trans- 
mission and reflection amplitudes. 



Transmission eigenvalues 

Our calculation yields 11 scattering channels (i.e. 11 
propagating states at the Fermi level) in the 5x5 atom 
Au(100) quantum wire, consistent with previous calcula- 
tions [9]. We obtain the transmission eigenvalues by di- 
agonalizing the 11x11 transmission block of the scatter- 
ing matrix connecting the left and right leads, sl R . The 
eigenvalues are shown in Fig. SI versus the APC gap size. 
One transmission eigenvalue dominates the others at all 
gap sizes, and becomes perfectly conducting when the 
atomic gap size is near the Au-Au interatomic distance 
of 2.87 A. We also find two sets of two-fold degenerate 
transmission eigenstates (red triangles) for all gap sizes. 



Damping and force noise 

To calculate the effective temperature plotted in Fig. 
2c of the main text, we first calculate the symmetrized 
force noise and damping from the scattering matrix using 
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FIG. SI: Transmission eigenvalues versus APC gap size. Red 
triangles indicates two-fold degenerate eigenvalues. The verti- 
cal dashed line marks the Au-Au interatomic distance of 2.87 
A in bulk gold. 




FIG. S2: Damping and force noise versus APC gap size, cal- 
culated from the scattering matrix obtained using DFT. Dips 
appear at the Au-Au interatomic distance of 2.87 A in bulk 
gold (vertical dashed line) in both quantities calculated from 
the ^-dependence of transmission only; however, no dips are 
present when phase corrections are included. Both 7 and Sf 
are scaled by their maximum values (from the full calcula- 
tion) . 



Eqs. (3) and (4) of the main text. T e g is then obtained 
from T e ff = Sf/ZM'j as discussed in the main text. The 
force noise and damping are shown in Fig. S2. For com- 
parison, we also show the force noise and damping cal- 
culated from the ^-dependence of the transmission only, 
neglecting the phase contributions. Both quantities cal- 
culated from T(x) alone show dips at a gap size near 
the Au-Au interatomic distance in bulk gold. This is not 
surprising, since this point corresponds to a minimum in 
the total free energy of the atomic system; moreover, the 
dominant transmission eigenvalue reaches a maximum at 
this point (see Fig. SI). Interestingly, corrections from 
the scattering phases completely wash out these dips. 
Note that the dips are not reflected in T e g; from the 
transmission terms only, we always obtain T e g = eV/2. 
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